Atom number filter in an optical lattice 
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We present an efficient procedure to filter out from an optical lattice, having inhomogeneous site 
occupation number, only preselected number of bosonic atoms per site and place them into another 
internal atomic state, creating thereby a lattice with desired site occupation number. 
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Ultracold atoms in optical lattices [T| represent a re- 
markably clean and controllable system [2, 3] to realize 
the fundamental Bose- Hubbard model [|[. Its two main 
ingredients are the atom tunneling, or hopping J, be- 
tween the neighbouring lattice sites and the on-site atom- 
atom interaction U. In a homogeneous lattice, when the 
kinetic energy due to the inter-site hopping dominates, 
J > U, the atoms are delocalized over the entire lattice 
yielding a superfluid (SF) phase, while in the opposite 
regime of strong on-site interaction, U 3> J, the hopping 
is energetically suppressed resulting in a Mott insulator 
(MI) phase with fixed integer number n of localized atoms 
at each lattice site. When a deep optical lattice is super- 
imposed by a shallow confining potential, there can be MI 
phases with occupation numbers of n = 0, 1, 2, ... in suc- 
cessive spatial shells [E|-0], separated by SF phases with 
intermediate mean occupation number corresponding to 
delocalized atoms on top of the filled MI shell. 

Experimentally Q, the quantum phase transition be- 
tween the SF and MI phases is implemented by adia- 
batically increasing the lattice depth which results in the 
reduction of intersite tunneling amplitude and simultane- 
ous increase of the on-site interaction |2J. If, however, the 
lattice potential is raised quickly, so that the tunneling 
is suddenly switched off, each site occupation "freezes" 
to whatever atom-number distribution it corresponded 
to just before the switching off, be it a SF, a MI, or a 
spatially-dependent combination of the two phases. 

In this paper, we propose a very efficient method to 
filter out from such a frozen ( J = 0) optical lattice only 
the desired number N of atoms per site. This is achieved 
by using an external field which couples the initially pop- 
ulated internal atomic state \a) to another internal state 
\b) trapped by a second optical lattice potential. We show 
that, for strong enough state- (or lattice-) dependent on- 
site interactions, the coupling field with properly tuned 
frequency will selectively transfer to the second lattice 
only the singles (N = 1), the pairs (N = 2), or the triples 
(N = 3) of atoms, via the corresponding TV-photon res- 
onant transition. Hence, after the transfer, the second 
lattice will only have the desired site occupation number 
N = 1,2, or 3, while the first lattice will contain all the 
other occupation numbers n N. 

Before proceeding, we note related, but different, ear- 
lier work. Rabl et al. @ proposed to reduce the site oc- 
cupation number defects in an optical lattice by adiabat- 
ically transferring a chosen number of atoms to another 



internal state. DeMarco et al. [5] studied similar sys- 
tems employing rapid adiabatic transfer of atoms to the 
second internal state, or inducing resonant single-photon 
Rabi oscillations between the atomic states with occupa- 
tion number-dependent Rabi frequencies. Mohring et al. 
@ discussed coherent extraction of atoms from a BEC 
reservoir into the quantum tweezers — tight trap — using 
adiabatic and resonant transfer techniques. 

Considering only two internal atomic states and cor- 
responding optical lattice potentials, the Hamiltonian of 
the system takes the form 

H = [{%u3 a + £a,j)h a ,j + ^U aa h a ,j(n a ,j - 1) 
j 

+ (huj b + eb,j)fib,j + \Ufohb,j{fib,j - 1) 
+U ab h a ,jh btj + mibfee-™ + &yb je iut )) . (1) 

Here a,j (aj) and bj (b^) are the annihilation (creation) 
operators for bosonic atoms in the internal states \a) and 
| b), of energies huj a and tkob, localized at lattice site j, 
with single-particle energies e a j and £bj, and fi a j = a^hj 

and hbj = bjbj are the corresponding number opera- 
tors. A natural basis for Hamiltonian (fTJ) is that of the 
eigenstates \n a ,j) of operators n a j whose eigenvalues 
n = 0, 1, 2, . . . denote the number of atoms in the cor- 
responding state | a) (a = a, b) at site j. Next, U aa — 
9aa J d 3 r\w a (r)\ is the on-site interaction energy for the 
atoms in state \a), and U ab = g ab J d 3 r|w a (r)| 2 |w b (r)| 2 
is the interaction between the \a) and \b) atoms, where 
g aa i = Ana aa ifi 2 /M , with a aa > being the corresponding 
s-wave scattering length, M the atomic mass and w a (r) 
the (localized) Wannier function of the lowest Bloch band 
of the corresponding lattice potential 0, 0. Finally, 
= Qab J d 3 rwl(r)wb(r) is the coupling amplitude be- 
tween the localized atoms in states \a) and |6), which is 
induced by an external field with the "bare" (free-atom) 
Rabi frequency £l a b- This field can be a microwave field 
of frequency w ~ u; b — uj a coupling the atomic hyperfine 
states | a) and \b) through a magnetic dipole transition, 
or an optical bi-chromatic field inducing Raman transi- 
tion |a) — > | b), in which case uj is the frequency difference 
between the two field components (the corresponding dif- 
ferential ac Stark shift of \a) and \b) can be incorporated 
in u> a or u>b). Note that the rotating- wave approximation, 
requiring -c w, is presumed in the last term of Eq. fl}. 
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In deep optical lattices, the Wannier functions w a (r — 
r a j) localized on individual sites j can be well approxi- 
mated [l| by the ground-state wavefunction of a harmonic 
oscillator centered at r Q 



3 > 



w a (r - r aJ ) 



3/4 



exp 



2u 2 



(2) 



where the width er Q = ^Jh/Mv a is expressed through the 
vibrational frequency v a = y / 2ir 2 V a /Md 2 determined by 
the lattice potential amplitude V a and period d. For the 
interaction parameters of Hamiltonian (JlJ we then obtain 
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where Sr = \r a j — Tf,,j\ < d is a possible offset of the 
lattice potentials for the atoms in states |o) and \b) 
[To| ■ These expressions attest to the controllability of 
the atom-atom interactions U aa > and coupling 51 through 
the interatomic scattering lengths a aa i (a, a 1 = a, 6); the 
optical lattice parameters, including the lattice modula- 
tion depths V a , affecting <t q 's, and the relative offset Sr; 
as well as the external coupling field amplitude, affecting 
Qah. In the experiments, typically U aa > /Ti < 2tt x 10 kHz 
U S B • We emphasise that the single-band approx- 
imation inherent in Hamiltonian (TIJ requires that the 
atom-atom interactions be small compared to the excited 
band energies, U aa > < Tw a . 

While, in general, e aj and Eb.j need not be uniform 
throughout the lattice, due to, e.g., shallow external trap, 
we assume that the difference \e a ,j — Sb,j\ is constant, and 
typically small compared to LJb — w a , for all j. Accord- 
ingly, we define TiuJba = huJb+Sb,j — (fu-o a +£a.j) = const V? 
and omit the subscript j from now on. 

Assume that initially all the atoms are in state a), 
with the sites of the corresponding lattice having arbi- 
trary occupations \n a ), n = 0, 1, 2, . . ., and all the sites 
of the other lattice empty, |0{,). 

Before describing our main idea of atom-number fil- 
ter, we briefly consider a simple but instructive case of 
uniformly interacting U aa i — U {a, a' — a,b), or non- 
interacting U = 0, atoms subject to a resonant cou- 
pling field ui = ojba- Within an iV-atom subspace (N = 
n a +nb), the transition matrix element of Hamiltonian (fTJ) 
between any pair of states of the form | (N — n) a , rib) and 
\{N -n - l) a , (n + l) fc ) is given by tt^f(N -n)(n+l), 
as dictated by the bosonic nature of the atoms. This 
coupling pattern makes the system formally analogous 
to a spin-J' in a magnetic field. Indeed, as we may 
recall from the theory of angular momentum, with the 
quantization direction along an axis perpendicular to 
the magnetic field direction, the matrix elements for the 



transitions 1^7,™} — > \ J : m+ 1} between the neighbor- 
ing magnetic sub-states (to = —J,...,J) are given by 
Sly/(J -m)(J + m + The spin then exhibits non- 
dispersive precession about the field direction with the 
Larmor frequency il. Setting formally J = iiV and 
to = n — leads to the above matrix elements of (JlJ. 
A curious consequence of this analogy is that, within any 
Af-atom subspace, the resonant coupling field will induce 
oscillations between states \N a ,0b) and \0 a ,Nb) with the 
same frequency J7. And in particular, starting from all 
the atoms in state \a) and arbitrary site occupation num- 
bers of the corresponding lattice, at time r = 7r/2£7, all 
the atoms will simultaneously be transferred to state |6). 

We now discuss the transfer of selected number of 
atoms N between the two lattices, as illustrated in Fig.Q] 
(left panel). This atom- number filtering procedure is very 
simple yet remarkably efficient and robust, provided 



\Uab- \{Uaa + Ubb)\ » «2, 

\u aa -u bb \, \u ab - u aa , bb \ > ha 



(4a) 
(4b) 



The first of these conditions ensures that within the se- 
lected A^-atom subspace all the intermediate states are 
nonresonant, while the remaining conditions are needed 
to suppress all the transitions out of the other initial 
states \n a ) with n ^ N, as clarified below. For conve- 
nience we denote SU = [U a b ~ \ {U aa + U b b)]/Tl. 

(i) Single atom transfer, N = 1. To filter out only 
the single atoms per site, we tune the frequency of the 
coupling field to be resonant with the atomic transition 
| a) -+ 1 6), i.e., we set u> = uj b a- The field will then induce 
resonant Rabi oscillations between the states |l a , 0;,) and 
|0 a , lb) with frequency fiW = £1. If we apply the field for 
time = 7r/2il^ 1 \ resulting in a 7r-pulse, all the single 
atoms |1 Q ) will be transferred to |1&). 

(ii) Two atom transfer, N — 2. To filter out only 
the pairs of atoms per site, we choose the frequency 
of the coupling field according to the condition 2lu = 
2ujba + (U b b — U aa )/h, which implies a two-atom (and 
two-photon) transition |2 a ,0fc) — > |l a ,l&) — > |0 o ,2{,} via 
nonresonant intermediate state |1 , 1&) detuned by SU. 
The corresponding two-atom (-photon) Rabi frequency is 
then f2(2) = 2V? jSU (the factor of 2 = 2! originates from 
double application of bosonic operators b^d to the initial 
state |2 a ,0h)), and at time r^ 2 ) = n/2Q( 2 \ correspond- 
ing to an effective 7r-pulse, all the pairs of atoms \2 a ) will 
be transferred to |2&). 

(Hi) Three atom transfer, N = 3. To filter out only 
the triples of atoms per site, we choose the frequency 
of the coupling field according to the condition 3oj — 
Suiba + 3(L/f,i> — U aa )/h, which implies a three-atom (- 
photon) transition |3 o ,0&) -> |2 OJ 1&) -> |l a ,2 h ) -> 
|0 Q ,3b) via nonresonant intermediate states |2 a ,lf,) and 
|l Q ,2b) both detuned by the equal amount of 2SU. The 
corresponding three-atom (-photon) Rabi frequency is 
then = 6fl 3 /(2SU) 2 (the factor of 6 = 3! originates 
from triple application of b'd to state |3 a ,0&)). Note 
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FIG. 1: Schematics (left column) and dynamics (right column) of atom transfer between two optical lattices: selective transfer 
of (i) single atoms, N = 1; (ii) pairs of atoms, N = 2; and (iii) triples of atoms, N = 3, via the corresponding iV-photon 
resonant transitions. The main graphs display the probabilities Pjv a for N atoms in the corresponding states \a) (a — a,b), 
while the insets show the probabilities P„ a of initial states \n a ) with n ^ N. (For n > 4, the probabilities P Ua oscillate with 
even smaller amplitudes and therefore not shown.) The numerical simulations employ the parameters U aa i and f2 listed in the 
text, the time is in units of f2 _1 and the evolution terminates at the corresponding t^ n '. 



that since both intermediate states |2 a , 1;,) and \l a ,2 b ) 
have the same detuning 2SU, the second-order ac Stark 
shifts of states 3 a ,0b) and a ,3f,} are the same, given 
by 3S1 2 / (2(5C7) , and the differential shift on the three- 
photon transition |3 a ,0&) — > |0 a ,3b) vanishes. Hence, 
applying the field for time r' 3 ' = tt/2Q( 3 \ corresponding 
to an effective 7r-pulse, all the triples of atoms |3 a ) will 
be transferred to |3f,}. 

The above procedure can be generalized to multipho- 
ton transfer of any number of atoms N between the 
two lattices. Under the iV-photon resonance condition 
Nlj = Nuj ba + ±N(N - l){U bb - U aa )/h, the effective 
TV-atom (-photon) Rabi frequency is then given by 



N 



[(N - 1)\} 2 SU N - 1 (N-l)\ 



n 

W 



N-l 



However, due to the above scaling of fiW and condi- 
tion f2 <§C \5U\, the corresponding transfer time = 
%/2Qy f > will become prohibitively long for N > 4 in a 
realistic optical lattice experiment, as discussed below. 

In Fig. [T] (right panel) we demonstrate, via numerical 
solution of the corresponding Schrodinger equations, that 
the transfer of the selected number of atoms N = 1, 2 and 
3 between the two lattices is indeed very efficient, with 
the probabilities iVjO 7 " ) °f the final states \Nb) at 
the corresponding times t^' being close to unity, while 
the probabilities P Ua of initial states \n a ) with n ^ N 
changing very little during the transfer. For these sim- 
ulations, we choose U aa /h ~ 2tt x 10 4 s _1 , and, upon 
assuming a aa — a>bb — o-abi V a /Vb = 3 and Sr — 0, obtain 
from Eqs. 43a} U bb ~ 0.44f7 aa and U ab ~ 0.63U aa . We 
then have 5U ~ 2tt x 930 s _1 and set fl = 2ir x 100 s _1 . 
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FIG. 2: Filtering out of the initial Poisson distribution with 
mean atom number (n) a = 1.5 (top panel graph) single atoms 
(i), pairs of atoms (ii), and triples of atoms (iii). In the lower 
panel graphs, the empty (red) bars show the atom-number 
probability distributions Pn a (T^ ^) remaining in the first lat- 
tice after the transfer. The transfer fidelities F(N) are defined 
in Eq. ©. All parameters are as in Fig. [T] 



The corresponding one-, two- and three- atom transfer 
times are given, respectively, by r' 1 ) = 2.38 x 10~ 3 s, 
t^ 2 ) = 1.16 x 10~ 2 s and = 0.144s, which are shorter 
than the typical lifetimes (0.5 s) of cold atoms in optical 
lattice MI shells with n<3|. 

Hence, using our procedure one can separate the spa- 
cial MI shells of the optical lattice with different atom 
numbers n @, @] , placing in another lattice only the de- 
sired A^th shell, which can be a filled sphere (or circle 



in 2D), or a hollow sphere (ring in 2D), depending on 
whether it is extracted from the central part of the trap 
or not. This is then followed by discarding (releasing) 
the atoms of the first lattice. Another useful application 
of our atom-number filtering technique is a preparation 
of pure samples of the interaction-bound lattice dimers 
[Ill4l4| . or trimers [15j . without resorting to more com- 
plicated procedures involving Feshbach association, pu- 
rification and dissociation of atom pairs [ll|, Gil • 

As an example, in Fig. [5] we illustrate the filtering out 
of the initial Poisson atom- number distribution P Ua (0) = 
(n)"e~( n '/n.!, corresponding to a frozen SF phase with 
mean occupation number {n) a = 1.5, the desired number 
of atoms N = 1, 2 or 3. The only variables adjusted to 
each N case are the coupling field frequency lj and trans- 
fer time N ^ , with all the other parameters the same, as 
described above in connection with Fig. [TJ We quantify 
the transfer using the fidelity 



F(N) 



E„>o^(rW)' 



(5) 



for which we obtain very high values F(l) — 0.999, 
F(2) = 0.987 and F(3) = 0.958. 

To conclude, we have proposed and analysed a very ef- 
ficient and robust procedure to filter out from an optical 
lattice with an arbitrary inhomogeneous site occupation 
number only preselected number of bosonic atoms per 
site and place them into another internal atomic state, 
creating thereby a lattice with desired site occupation 
number, which we envision to have a number of interest- 
ing applications. 
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